ON THE PROPAGATION OF OCEANIC WAVES DRIVEN BY A 
STRONG MACROSCOPIC FLOW 



ISABELLE GALLAGHER, THIERRY PAUL, AND LAURE SAINT-RAYMOND 



Abstract. In this work we study oceanic waves in a shallow water flow subject to strong 
wind forcing and rotation, and linearized around a inhomogeneous (non zonal) stationary 
profile. This extends the study Q], where the profile was assumed to be zonal only and where 
explicit calculations were made possible due to the ID setting. 

Here the diagonalization of the system, which allows to identify Rossby and Poincare 
waves, is proved by an abstract semi-classical approach. The dispersion of Poincare waves 
is also obtained by a more abstract and more robust method using Mourre estimates. Only 
some partial results however are obtained concerning the Rossby propagation, as the two 
dimensional setting complicates very much the study of the dynamical system. 



1. Introduction 

This paper is a continuation of pQ so before discussing the matter of this paper (in Section [2]) 
let us review the contents of that work. We shall start by recalling briefly the model, then we 
shall explain the methods and results obtained in [T] and discuss their limitations. 

1.1. The model. The goal of [lj is to understand, through the study of a toy model, the 
persistence of oceanic eddies observed long past by physicists among which [H El El El El ED] , 
who gave heuristic arguments to explain their formation due both to wind forcing and to 
convection by a macroscopic current. 

The ocean is considered in this toy model as an incompressible, inviscid fluid with free sur- 
face submitted to gravitation and wind forcing, and we further make the following classical 
assumptions: we assume that the density of the fluid is homogeneous p = po = constant, that 
the pressure law is given by the hydrostatic approximation p = pogz, and that the motion 
is essentially horizontal and does not depend on the vertical coordinate. This leads to the 
so-called shallow water approximation. 

For the sake of simplicity, the effects of the interaction with the boundaries are not discussed 
and the model is purely horizontal with the longitude x\ and the latitude X2 both in R. 

The evolution of the water height h and velocity v is then governed by the shallow-water 
equations with Coriolis force 

d t { Po h) + V -(pohv) = 

(1.1) 

dt(pohv) + V ■ (pohv ®v) + uj(pohv) + p ghVh = p hr 
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where w denotes the vertical component of the Earth rotation vector Q, v 1 - := (—V2,vi), g 
is the gravity and r is the - stationary - forcing responsible for the macroscopic flow. The 
vertical component of the Earth rotation is therefore Vt sin(x2/R), where R is the radius of 
the Earth; note that it is classical in the physical literature to consider the linearization of u 
(known as the betaplane approximation) oo(x2) = £1x2/ R- We consider general functions ui in 
the sequel, with some restrictions that are be made precise later. 

We consider small fluctuations around the stationary solution (h, v) satisfying 

h = constant, V • (v (g) v) + ujv 1 - = r, div v = 0. 

In pQ the study is restricted to the case of a shear flow, in the sense that v(x) = {vi{x2)-, 0), 
with vi a smooth, compactly supported function. Moreover some orders of magnitude and 
scalings allow to transform the previous system into the following one: 



(1.2) 



dtr\ H — V ■ u + u ■ 



dtu + -^bu 1 - + -Vr? + u ■ Vu + u ■ Vu + e 2 u ■ Vu = , 

e z e 



where b := and where e is a small parameter (of the order of Fr 2 , where Fr is the Froude 

number, and of R02 where Ro is the Rossby number). 

1.2. Methods and results in pQ. Most of the analysis in [l] concerns the linear version 
of (jl.2p . namely the following system: 

(1.3) e 2 id t v + A(x2,eD,e)v = v = (v , v\, V2), 

where D := id, and the linear propagator is given by 

(eu\edi ed\ £$2 
ed\ eu\ed\ —b{x2) + £ 2 u[ 
E&2 b(x2) euiedi 

The first step of the analysis consists in diagonalizing (approximately) the system (jl.3p . The 
computation of a kind of characteristic polynomial associated with (|1 .3[) . in symbolic form, 
allows to construct three symbols the quantization of which provides three scalar propagators 
(this will be explained more explicitly below). 

Two of those propagators, called Poincare propagators, are then proved to satisfy dispersive 
estimates; that result relies on a spectral analysis (usual semi-classical theory does not operate 
here due to the very large time scales at play) using Bohr-Sommerfeld quantization, which 
requires that b has only one, non degenerate critical value and which also uses very much the 
fact that the motion is translation-invariant in x\. A stationary phase argument on the 
spectral decomposition of any solution to the Poincare propagation gives the result: Poincare 
modes exit any compact set in finite time. 

The last propagator is the Rossby one, which is one order of magnitude (in e) smaller than 
the Poincare modes. This allows to analyse the propagation by semi-classical analysis tools. 
In particular the precise study of the dynamical system associated with those waves, which is 
an integrable system due to translation invariance in x±, allows to derive a condition on the 
initial microlocalization of the solution which guarantees that the Rossby waves are trapped 
for all times in a compact set. 
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Those results on the linear system (II. 3p can finally be transposed to the original system (|1.2p 
due to the high power of e in front of the nonlinearity, and due to the semi-classical set- 
ting, which allows to exhibit vector fields which almost-commute with the linear opera- 
tor A(x2, sD, e). 

1.3. Limitations of the methods of [I]. The restriction which is the most used in the 
analysis described briefly in the previous paragraph is the fact that the stationary flow u is a 
shear flow of the type u = (ui(x2), 0). Indeed 

• It allows to Fourier-transform in the direction x%, which makes the diagonalization 
procedure much easier; 

• It simplifies the spectral analysis of Poincare waves, again due to the Fourier transform 
(in particular the dual variable £i is fixed during the propagation, and there is a wave- 
like behaviour in xi); 

• It allows the Rossby dynamical system to be integrable, which is a tremendous help 
in the analysis. 

An additional restriction in the previous arguments is that in order to prove the dispersion of 
Poincare waves, the rotation amplitude b should have only one, non degenerate critical value: 
this allows to use a Bohr-Sommerfeld quantization argument to compute the eigenvalues of 
the Poincare operator. This assumption on b is not really restrictive from the physical point 
of view. On the other hand, it is important for physical reasons to consider 2D convection 
flows. 

1.4. On the nonlinear term. As explained above, most of the analysis in [T] is concerned 
with the linear system (jl.3j) . In order to transpose the linear results to the nonlinear setting, 
one uses the following arguments (along with the fact that the coupling is vanishes when e 
goes to zero): 

• Uniform existence which is obtained via an almost-commutation result; 

• Bilinear estimates in anisotropic semi-classical spaces; 

• A Gronwall lemma, which requires an L°°(R 2 ) bound on the linear solution. This is 
not known in general, due to the bad Sobolev embeddings in semi-classical settings, 
so the nonlinear result is proved for vanishing couplings only. 

It is important to notice that none of those three steps require that u is a shear flow. In 
the whole of this paper we shall therefore only focus on the linear equation, and leave to the 
reader the transposition to the nonlinear equation, using the above steps. 



2. Main result of this paper and strategy of the proof 
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We shall assume throughout the paper that b is smooth, with a symbol-like behaviour: for 
all a € N, there is a constant C a such that for all x 2 G R, 

(2.3) \b^(x2)\<C a {l + b 2 (x 2 ))K 
We shall further assume that 

lim b 2 (x2) = oo, 

|x 2 |— >oo 

and that b 2 has only non degenerate critical points. 

We shall also suppose that the initial data is microlocalized in some compact set C of T*R 2 
satisfying 

(2.4) C n{e i +e 2 +b 2 (x 2 ) = 0} = ® 
or actually rather 

(2.5) Cn{£i = O} = 0. 

We shall prove that assumption (|2.4p is propagated by the flow, while ()2.5[) is propagated only 
by the Poincare component. We recall (see for instance [I], Appendix B) that a function / is 
microlocalized in a compact set C of T*R 2 if for any (xo,£o) in the complement of C in R 4 
(we shall identify T*R 2 to R in the following), there is a smooth function xo> bounded as 
well as all its derivatives and equal to one at (xo,£o)> satisfying 

(2-6) ||Opf(xoHllL 2 (R 2 ) = 0(^ 00 ), 

where Op^ 7 denotes the Weyl quantization: 



(2.7) Opf (xoM*) := 7^)4 y e i(a; - s/K/e Xo(^,0^(y) «. 

We also recall that (12 . 6[) means that for any N € N, there are eo and C such that 

Ve e]0,e o ], l|Op e H/ (xo)«o||L2(R 2 ) < C^. 

In the following, to simplify some formulations, we shall denote by (/^Supp*/ the projection 
of the (micro)support of / onto the * = axis. 



2.2. Statement of the main result and organization of the paper. Let us state the 
main theorem proved in this paper. 

Theorem 1. Let v £i o be a family of initial data, microlocalized in a compact set C satisfy- 
ing Assumptions \2.J$ - $2^fy - For any parameter e > 0, denote by v e the associate solution 
to \2. 1}) . Then for all t > one can write v e (i) as the sum of a "Rossby" vector field and a 
"Poincare" vector field: v £ (i) = v^(t) +v£(t), satisfying the following properties: 

(1) nSuppvf(t) and ^,Suppvf(t) satisfy \2.J$ for all times. 

(2) For any compact set 17 in R 2 , one has 

Vt>0, ||vf(t)|| L2(n) =0( e 00 ). 

(3) jiSupp X2 vf{t) lies in a bounded subset o/R uniformly in time. 
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Compared to [I], the main difficulties are due to the presence of a xi-dependent underly- 
ing flow u. The diagonalization of the system (exhibiting Rossby and Poincare-type waves, 
with very different qualitative features) must be revised, and obtained in a less explicit way. 
Moreover the proof of ([2]) in Theorem [H namely the dispersion of Poincare waves can also 
not be proved in the same way (note that it is not assumed here that b 2 has only one non 
degenerate critical value) . Finally the trapping of Rossby waves seems much harder to obtain 
since the underlying dynamical system no more decouples; the behaviour of the Rossby waves 
is therefore much less precise than in pQ. 

Let us explain our strategy here, compared with that in [1] described above. 

2.2.1. The diagonalization. The construction of the Rossby and Poincare modes is not as di- 
rect as in pQ due to the lack of translation invariance in x\. We choose therefore to follow 
a more abstract way to recover those modes in Section [3l which relies on semi-classical 
analysis, and normal forms (instead of explicit computations as in pQ). Finding the prop- 
agators associated with those modes requires a microlocalisation assumption of the type (|2.4p . 
in order for the eigenvalues of the matrix of principal symbols to be well separated. The diag- 
onalization result is therefore in this paragraph conditional to the fact that the solution to the 
propagation equation is correctly microlocalized (that corresponds to Point Q] of Theorem [T]) . 

2.2.2. Dispersion of Poincare waves and propagation of the nondegeneracy assumption \2. 5\) . 
In order to prove ([2]) in Theorem [1] we again rely on a more abstract, and more efficient 
method than that followed in pQ. It is based on Mourre estimates and the assumption (|2.5p 
on the initial data: we start by proving, by a semi-classical argument, that after a very short 
time (of the order of e) the support in x\ of the solution escapes the support of u. Then 
we use Mourre estimates to prove that the solution remains outside the support of u for 
all times, and actually escapes any compact set in x\ in finite time (to prove this last point 
we use the fact that the equation reduces to a translation-invariant equation in x\ since the 
support of the solution is outside the support in x\ of u). This allows finally to check that the 
nondegeneracy assumption (12. 5p does hold for all times. This analysis is achieved in Section HI 

2.2.3. Study of Rossby waves and propagation of the nondegeneracy assumption \2.J$ . In Sec- 
tion[5]we first prove that the nondegeneracy assumption (|2.4p does hold during the propagation 
of Rossby waves. That is due to semi-classical analysis, by the study of the dynamical system 
associated with those waves. The study of that system is also the key to the proof of Point (|3|), 
which is also proved in Section [5j 

3. Reduction to scalar propagators 

In this section we shall construct three operators T+, T_ and Tr diagonalizing A(x,eD,e). 

We shall start by proving a general diagonalization result, and at the end we shall apply the 
general result to our context. 

Before stating the general result, let us give some notation. A semi-classical symbol is a 
function a = a(x,£;e) defined on R M x]0,eo] for some e$ > 0, which depends smoothly 
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on and such that for any a £ N M and any compact /C C R 2c( , there is a constant C 

such that for any 6 /Cx]0, £o], 

ia a a((^e)^)l <c. 

We shall consider the Weyl quantization of such symbols, as recalled in (|2.7p : for all u is 
in V(R d ), 

Opf (a)u(x) := | e^^M^, «. 

We shall say that a pseudodifferential operator Op^(a) is supported in a set /C if for any 
smooth function x equal to one in a neighborhood of /C one has ax = a. 

Finally we shall say that a matrix is pseudodifferential if each of its entries is a pseudodiffer- 
ential operator. 

Let us first prove the following general result. 

Theorem 2. Let K, be a compact subset of~R 2d , and consider a N x N hermitian pseudodif- 
ferential matrix A e = A(x,eD,e), supported in K. Assume that 

• the (matrix) principal symbol of A(x, eD, 0), denoted by Ao, is diagonalizable, in the 
sense that there are some unitary and diagonal matrices of symbols, U and T>, such 
that 

U^AoU = V, 

• the eigenvalues (Si(x, £), . . . , 6n(x, £)) satisfy 

(3.1) Vi^i, inf \S i {x^)-S j (x,0\>C>0. 

(o:,£)e/C 

Then there exists a family of unitary and diagonal pseudodifferential operators V e and D £ 
supported in K, such that: 

(3.2) V*A £ V e = D £ + 0(e°°), V*V £ = 1 + 0(e°°). 
Moreover one has 

(3.3) D £ = Opf (X>) + eDx + 0(e 2 ), 
where the principal symbol of D\ is given by 



T>\ = diag f Ai - ^ IlD ° ^j 



with the notations 



A a = - (OpfCZTKOpf (W) - D ) , 
(3.4) e 

/^-(Opf^Opf^)-/) 

More explicitly, let us denote by aij(x,^) the matrix elements of A\, subsymbol of A{x,eD,e) 
defined by: 

OA 

At(x,0 i=-£(x,£,0) 
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and by u n j(x,£), i = 1 . . . N, the coordinates of any unit eigenvector of Ao(x, £) of eigenvalue 
5 n (x,£). We have 

(3.5) (V 1 ) nn = (%(u^{a jk ,u kn }) + My*n> \ + {WAlU ) nn , 

jk=l...N ^ ' ' 

where {f,g} ■= V^/V^g — V^/V^g «s the Poisson bracket on T*R". 

Here and in all the sequel, we say that a pseudo-differential operator V is unitary if it satisfies 

V*V = 1 + 0(6°°). 

The proof is divided into two parts: in Section [3.11 we present the formal construction and in 
Section [JO] we show that the symbols of the various operators formally constructed are indeed 
symbols. Finally Section [3.31 is devoted to the case of the matrix given by (|2.2p . 

3.1. The formal construction. The proof of Theorem [2] is a combination of semiclassical 
and perturbation methods. Let us start by defining 

U = OpY(U). 

Elementary properties of the Weyl quantization imply then that UqA £ Uo = Dq + 0(e). 

The following proposition shows that one can construct a unitary pseudo differential opera- 
tor Uoo such that 

U^A £ U oo = D + O(e). 
Lemma 3.1. Let U be a pseudodifferential matrix such that U*U = I + el\, where I is the 

oo 

identity. Then one can find V ~ ^^ £fc ^4 such that 

k=0 

(3.6) (U + eV)*(U + eV) = I + 0(e°°). 

Proof. Let us denote V := ~^hU. On easily check that (U + eV )*(U + eV ) = I + 0(e 2 ). 
Indeed 

(U + sV )*(U + eV ) = U*U + ^(hU*U + U*Uh) + 0(e 2 ) 

= l + eh-eh + 0(e 2 ). 
Then one concludes by iteration. □ 

That lemma allows to define the pseudo-differential operator of (semiclassical) order 

Ai = i (U^A e Uoo - D ) , 

where Uoo is a unitary operator. 

Now our aim is to find a unitary operator Voo (up to 0(e°°)) such that 

(UooVooYA^UooVoo) = Doo + 0(e°°), 
where = Dq + eD\ + . . . is a diagonal matrix satisfying the conclusions of the theorem. 
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We shall write = e teW , with W selfadjoint (so thus constructed is automatically 
unitary). We recall that if W is a pseudodifferential operator, then so is e ieW (simply by 

writing e ieW ~ J2o ^W k ). 

oo 

We look for W under the form W ~ ' s ^e k Wk 1 and compute the W k recursively. Since 

o 

(is) 2 

V^(D + eAijV^ = (Do + eAx) + ie[(A, + eAi), W] + ^H[(D + eAx), W], W] + . . . 
we see that, if W\ satisfies 

(3.7) [Do, Wx] + Ai = Di + 0(e 2 ), D x diagonal, 
then we have that 

(3.8) e- ieW '(D + eAx)e ieW ' =D + eD 1 + e 2 A 2 , 

where A 2 is a zero order pseudodifferential operator. The following lemma is a typical normal 
form type result, and is crucial for the following. 

Lemma 3.2. Let Dq be a diagonal pseudodifferential matrix whose principal symbol T>q has 
a spectrum satisfying \3. 1\) and let Ai be a pseudodifferential matrix. 



Then there exist two pseudodifferential matrices W and D\ , with D\ diagonal, such that: 
(3.9) [D ,W]+Ax = Dx + eA 2 , 

where A 2 is a pseudodifferential matrix of order 0. 

Moreover the principal symbol of D\ is the diagonal part of the principal symbol of A]_. 

Proof. By the non degeneracy condition of the spectrum of T>q we know, by standard ar- 
guments (see [11] for instance), that there exists a matrix Wo and a diagonal one T>\ such 
that 

[V ,Wo}+T>x,a = V 1 , 
where T>\ q is the principal symbol of Ai. 



Indeed it is enough to take T>\ as the diagonal part of T>\ § and 

(3.10) (m(x,o)i 



_ (T>x,o(x,Z))i,j 



and notice that the Weyl quantization of Wo satisfies (j3.9|) . □ 

By Lemma 13.21 we know that there exists Wx satisfying (13.T|) . Developing (13.8P as 

e~ ieW ^{D + eA 1 )e l£m = D + e(Ax + [£>„, Wx}) + e 2 A 2 , 
we get immediately (|3.8p . 

It is easy to get convinced that all the Wk will satisfy recursively an equation of the form 

[D ,W k ] + A k = D k + O(s), 

which can be solved by Lemma [ 
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The expression for the principal symbol of D\ follows by construction and the following well 
known lemma (see [8] for instance): 

Lemma 3.3. Let a and b two symbols. Then the principal symbol of Op^(a)Op]^(&) is ab 
and its subprincipal symbol is ^{a,b}. 

In order to derive (|3.5p we have to compute the subprincipal symbol of the diagonal part of 
the right-hand side of (|3.4p . that is, for each n = 1 . . . N, 

^Opf (Z7~)Opf ((.4o + sA^OpfiUkn), 

jk 

since U is unitary. 

The term eA\ is obviously responsible for the second term in the right-hand side of (|3.5p . 
Using Lemma [3 .31 and the distributivity of the Poisson bracket, we get the following expression 
for the first one: 

jk 

= ^ (Uj n {(Ao)jkMkn} + {At)) jk{UjnMkn} + ^kn{^jn, {-A))jk}) 
jk 1 

Interverting j and k in half of the terms and noticing that, since .4,0 is Hermitian, (Ao)jk = 
(<Ao)kj> we g et easily (pU^ . 

3.2. Symbolic properties. 

With the hypothesis that both T> and IA are pseudo differential matrices it is quite obvious 
that V £ and D £ are pseudodifferential matrices as well. Indeed the formal construction in the 
preceding section shows that the iterative process uses only three things: multiplications of 
pseudodifferential operators, computation of subprincipal symbols and solving equation (|3.9f) . 

For (|3,9p . the formula (|3.10p used in the proof of Lemma l3.2l together with the non-degeneracy 
condition (|3.ip which shows clearly that (<5j(x,£) — 6j(x, is a symbol, implies that Wo is 
a pseudodifferential operator. 

Note that the microlocalization assumption is crucial in order that the expansions obtained 
by this iterative construction do define symbols. We have indeed no uniform control on the 
growth at infinity. 

3.3. The Rossby-Poincare case. 

In the case of oceanic waves A(x,sD,s) is given by (|2.2p : 



A[x,eD) 




eu ■ eV + e 2 d\ui 
b(x 2 ) + e 2 diu 2 



ed 2 

-b(x 2 ) + e 2 d 2 ui 
eu ■ sV + £ 2 d 2 u 2 
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Therefore 

[Q £1 

and 

'u ■ i 

u-i 

u-ij 

A straightforward computation shows that the spectrum of Aq is 





h- > 







) =«•«! 





1 


• 




Vo 








{o, ^ 2 + e 2 2 + & 2 M, 



3.3.1. Microlocalization. The three eigenvalues of Ao are separated if and only if 

Therefore, considering a compact subset /C of R M such that 

/c n {(xi, x 2 , 6, 6) /e? + £l + 6 2 (x 2 ) = 0} = 

ensures that 

• the eigenvalues do not cross, so that it is possible to get a unitary diagonalizing matrix 
with regular entries; 

• the non degeneracy condition ()3. 1|) is satisfied. 

In other words, A(x, eD, e) satisfies the assumptions of Theorem [2] provided that one considers 
only its action on vector fields which are suitably microlocalized. 

We assume of course that this microlocalization condition is satisfied by the initial datum, 
which is the condition (12.411. 



Furthermore, we will prove in the next two sections that the propagation by the scalar oper- 
ators T± and Tr (to be defined now) preserves this suitable microlocalization, thus justifying 
a posteriori the diagonalization procedure for all times. 

3.3.2. Computation of the Poincare and Rossby Hamiltonians. The above computations show 
that one can define the two Poincare Hamiltonians as follows: 



and we shall denote the associate operator constructed via Theorem [2] by T±. 

Now let us consider the Rossby Hamiltonian. In all this paragraph, for the sake of readability, 
we will denote 



An easy computation shows that a (normalized) eigenvector of Aq(x,£) of zero eigenvalue is 
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By Theorem O the Rossby Hamiltonian is then given by the formula 

(3.11) t r = ^2 (%{u]5{a jk ,u k0 }) + a ^ u ^ Uk ^ \ + ^ {Ai) jk Ujou ko . 

j,k=l...3 ^ ' j,k=l...3 

In order to compute the different Lie brackets, we start with a couple of simple remarks : 

{^,f} = d Xj f and {b(x 2 ),f} = -ft'(x 2 )%/. 
In particular, if / does not depend on x%, then {gi, /} = 0. 

The contribution of the first term in the parenthesis in (|3,lip is 



(u j0 {a jk ,Uko}) 

j,k=l...3 



(0 b { 6) 


(0»i 




-ib£i d 

(0 b X2 


1 


IS 








(01 







(Ob J I ' (0 



J x 2 jp- + TTT - ^ TTT 



Using the distributivity of the Poisson brackets, we get the contribution of the second term 
in a very similar way 

E ajkjujo, Uko} 
2 

i,fc=1...3 

b i& 1 f ft *Cl 1 f *£l ^2 



(0& '(£>*, J I (Ob' (Oft J l(C) 6 '(0 



6 ( (I { + § { &) ^} + - H { 6) 4} - 1 { 



(02 

The computation of the second term of the right hand side of (13. lip is trivial since A\ is a 
multiple of the identity. Adding with the two previous expressions we get finally 

gift' , 

TR = e 1 + tl + K X2 y +u ^ 

and the associate operator will be denoted by Tr. 

Remark 3.4. Since the elementary steps of the diagonalization process use only multipli- 
cations, computations of subprincipal symbols and solving normal forms equations, all the 
subsymbols of Tr and T± depend on X\ only through u and its derivatives. 
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4. Study of the Poincare waves 
In Section [3] we constructed two linear operators, called T±, whose principal symbols are 

We now want to study the propagation equation associated to those operators, namely the 
linear equation in R x R 2 

(4.1) ie 2 dw± = T±<p±, tp±\ t=0 = <p°± 

where tp± are microlocalized in a compact set C satisfying Assumption (12. 4j) . Before studying 
that equation we need to check that it makes sense, since a priori T± is only defined on vector 
fields microlocalized on such a compact set. This is achieved in the coming section, where we 
check that the separation of eigenvalues (|3.ip required in the statement of Theorem [2] holds 

because v/Cf"+£| + b 2 {x2) remains bounded away from zero during the propagation. 

Then we shall show that the solutions to these equations exit any compact set in finite time 
(Point © of Theorem d]). 

4.1. Microlocalization. Let us prove the following result, which provides the first part of 
Point ([TJ in Theorem Q] and allows to make sense of Equation (|4,ip for all times. 

Proposition 4.1. Under the assumptions of Theorem^ the operators T± are self-adjoint, 
and the function <p(t) = e % ~^ T± Lp\ are such that /j,Suppip±(t) satisfies \2.J$ for all times. 

Proof. The proof of that result relies on a spectral argument. Due to the form of the principal 
symbols of T± recalled above, the operators T± are self-adjoint. We can therefore define two 
families (V'n Wn of (pseudo)-eigenvectors of T± in L 2 (R 2 ) and two sequences of eigenvalues A J 
such that if the initial data writes 

^(x) = ^c±'V±(x), 

n 

then 

<P±(t,x) =^ e '^ c nVn( a: )- 
n 

Since the eigenfunctions ijj^ are microlocalized on the energy surfaces of the Poincare Hamil- 
tonians, the result follows. □ 

4.2. Dispersion. In this paragraph we shall prove Point ([2]) of Theorem [TJ The strategy is 
the following. 

In Section ^. 2. ll we prove using semi-classical analysis that for a very short time, the solutions 
to (14. ip remain microlocalized in a compact set satisfying assumption (|2.5p . and such that 
^tSupp a . 1 (p± become disjoint from Supp^u. Section T4.2.2I is then devoted to the long-time 
behaviour of the solution, and Mourre estimates allow to prove that the solution exits any 
compact set after some time, and that it remains microlocalized far from £i = 0. 
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4.2.1. Short time behaviour. The aim of this paragraph is to prove the following result. It 
shows that the solutions of (|4.ip exit the support of u after a time i ex it£, for |t ex it| large 
enough (independent of e). We only state the forward in time result: the backwards result is 
identical, up to changing the sign of time. We shall further restrict the analysis to T + since 
the argument for T_ is identical, up to some sign changes. 

Proposition 4.2. Let (p° be a function, microlocalized in a compact set C satisfying As- 
sumption ^2.4^ , and let ip be the associate solution of fl^. Ity . Let [u_,u+] be a closed interval 
of R containing Supp Xl u. There exists a constant t ex it > such that for any e g]0, 1[, the 
function (p(et ex it,-) is microlocalized in a compact set IC such that the projection of K, onto 
the x\-axis does not interesect [u-,u + ~\. Moreover fxSupp^cp is unchanged. 

More precisely, if [iSupp^ipa C R + \{0}, then fj,Supp xl ip(st ex it, •) c]u+,+oo[, and if fiSupp^ipo C 
R~ \ {0}, then /j,Supp xl ip(et ex it, ■) C] — oo, u_[. 



Proof. Define the function t/j(s) := (p(es). Then (14. ip reads 
(4.2) ied,il> = r+^ ; V|,=o = <P°, 

and any result proved on i/j on [0, T] will yield the same result for ip on [0, Te]. Notice 
that (|4.2p is written in a semi-classical setting, so by the propagation of the microsupport 
theorem (see for instance [8] and the references therein), the microsupport of ip is propagated 
by the bicharacteristics, which are the integral curves of the principal symbol. Recall that the 
principal symbol of T + is 



T + (a 1 ,x 2 ,a 2 ) = y/e 1 +e 2 + b 2 (x 2 ) 

and the bicharacteristics are given by the following set of ODEs: 

i' = V € T + (^,4,^), X° = (X?,X°) 

e* = -v,r + (^,*u*), e = 

Notice that r + is independent of x±, so is identically zero and therefore £\ = £?■ So for 
all s > 0, the microlocal support in of ip(s) remains unchanged, and in particular is far 
from £i = 0. Moreover one has 

x\~ ^ 



V($) 2 +m 2 +v(xi) 

Now we recall that the bicharacteristic curves lie on energy surfaces, meaning that on each 
bicharacteristic, t + (^,x 2 ,C 2 ) is a constant. That implies that (£, 2 ) 2 + 6 2 (x^) is a constant on 
each bicharacteristic, so that for all times, 

x\ = 4l 



If > 0, then x\ is propagated to the right and eventually escapes to the right of the support 
in x\ of u, whereas if £® < 0, the converse (to the left) occurs. Proposition 14.21 is proved. □ 
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4.2.2. Long time behaviour. The aim of this paragraph is to prove the following result, which 
again is only proved for positive times for simplicity. 

Proposition 4.3. Under the assumptions of Proposition {4~!\ ^ be the solution of fl^. i| j 
associated with the data </?(etexit> •)• Then fj,Supp Xl tp + (t) does not intersect Supp Xi u for t > 
etex.it, and fj,Supp^ 1 ip + (t) remains unchanged for t > ei e xit- Finally fiSupp Xl <p + (t) exits any 
compact set in x\ in finite time. 

Proof. Before going into the proof, we shall simplify the analysis by only studying the case 
of T + (the case T_ is obtained by identical arguments), and we shall only deal with the case 
when the support in £i of the data lies in the positive half space. The other case is obtained 
similarly. 

The proof is based on Mourre's theory which we shall now briefly recall, and we refer to [3] 
and [2] for all details. Let us consider two self-adjoint operators H and A on a, Hilbert space T~L. 
We make the following assumptions: 

(1) the intersection of the domains of A and H is dense in the domain T>(H) of H . 

(2) t i-4 e itA maps V(H) to itself, and for all <p° G V(H), 

sup \\He itA f°\\ < oo. 
te[o,i] 

(3) The operator i[H, A] is bounded from below and closable, and the domain T>{Bi) 
where iB\ is its closure, contains T>{H). More generally for all n G N the opera- 
tor i[iB n , A] is bounded from below and closable and the domain T>{B n+ \) of its clo- 
sure iB n+ \ contains V(H), and finally B n+ \ extends to a bounded operator from T>(H) 
to its dual. 

(4) There exists 9 > and an open interval A of R such that if E A is the corresponding 
spectral projection of H, then 

(4.3) E A i[H,A]E A >9E A . 

Note that Assumptions ([T]-[3]) can be replaced by the fact that [f(H),A] is bounded for any 
smooth, compactly supported function / (see [2]). 

Under those assumptions, for any integer m G N and for any 9' g]0, 9[, there is a constant C 
such that 

|| X _( A -a- 9't)e~* m g(H) X+ (A - a)\\ < Ct~ m 

where \± ls the characteristic function of R 1 * 1 , g is any smooth compactly supported function 
in A, and the above bound is uniform in a G R. As pointed out in [2], this implies in particular 
that for any ip° in the image of E A , the function e~ lHt tp° has spectral support in [a + 19', oo[ 
with respect to A, up to t~°° . 

Let us apply this theory to our situation. We consider equation (|4.ip with data e*~ fi e it ' r +<p9., 
and let us define the operator T9 as the operator T + where u has been chosen identically zero. 
We shall start by studying the equation 

(4-4) ie 2 d t £ = T% ^ =rfcxit = e*^^, 

for which we shall prove Proposition 14.31 Then we shall prove that the solution </> actually 
solves the original equation f)4. 1[) with the same data e l ~ E ? LT+ tp\ at t = et ex it up to 0(e°°), 
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because its support in x\ lies outside the support of u and because the symbolic expansion 
of T+ depends on x\ only through u and its derivatives (see Remark 13. 4ft . 

So let us start by applying Mourre's theory to (|4.4p . Let us write the projection of K, onto 
the £i-axis as included in [c?o,di] with < do < d\ < oo. We recall that on the support 

of e 1 c e' tT+ tp+, x\ remains to the right of the support of u. Then we apply the theory to H = 7+ 
and to A = x\ (the pointwise multiplication). Assumptions (H|) to ([3]) are easy to check, in 
particular because this is a semiclassical setting, so only the principal symbols need to be 
computed. Similarly finding a lower bound for E A i[T+, x\\E A boils down to computing the 
Poisson bracket {r + ,xi} where 

{f,g} = Vef-Vz9-Vxf-Vt9, 

and one finds 

(4.5) {r+, Xl } L < 



Since T? has constant coefficients in x\, one can take the Fourier transform of (|4.4j) and £i is 
preserved, so in particular for all times one has /^Supp^ 1 ^(t) C [c?o,di]. One can furthermore 
choose for A an interval of R of the type ]Dq, Di[ where the constants Do and D\ are chosen 
so that for any (x,£) £ IC, one has 



(4.6) D <^e i +e 2 +b 2 {x 2 )<D 1 . 

As the microlocal supports of the eigenfunctions of T? lie on energy surfaces, we know that 
the solution to (|4.4j) will remain in E A for all times. Now let us apply the results of [3] and [2]. 
By Lemma [3T3l (14. 5j> . f)4.6|) and the assumption on £j written above, we have that 

E A i[H,A]E A >e^-E A , 

so ((33J) holds with 9 = ed /D . It follows that the solution e ^ + (e l ~ 1 + ip%) to (Oil 
has a support in x\ such that 

d t 

Xi > u + + — — 
D e 



which proves the result for (j4.4j) . 

Since /xSupp Xl e* ? T+ (e i_£ t liT ' + (/9 C {.) does not cross Supp^u, one has actually 

e J 2 + (e e y+J = e ^ ( e e f+) m ^ 
locally uniformly in t (see Appendix). 

The proposition follows. □ 



5. Propagation of the Rossby waves 



5.1. Semiclassical transport equations and microlocalization. 
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Because of the scaling of the Rossby hamiltonian (which is smaller than the Poincare hamil- 
tonians by one order of magnitude), on the times scales considered here, the propagation of 
energy by Rossby waves is described by the hamiltonian dynamics 

dxj _ 8tr d& __ dr R 
dt d^i ' dt dxi ' 

which can be written explicitly 

dX2 nl _/, 



(5.1) 

d£i 
dt 



2b'{x 2 )^ + u 2 {x) 
dt (Ot 



<9iui(x)£i - diu 2 (x)&, 
dt 2= 2b(b'Y- b"(t) l 

dt Kl d)i 



— = fi TTVi d 2Ul {x)^ - d 2 u 2 {x)t > 2 



where we recall that = + £| + b 2 (x 2 ) . In order for the dynamics to be well defined 
and also in order to justify the diagonalization process, we need the quantity to remain 
bounded from below for all times. 

Proposition 5.1. Let C be some compact subset o/R 4 such that 

Cn {(xi,x2,£i,6) /£? + f 2 + b 2 (x 2 ) = 0} = 0. 

Then the bicharacteristics t \— )■ (x(t),£(t)) of the Rossby Hamiltonian starting from any point 
Ci > £2) of C are defined globally in time, and Vi G R, 

inf (£ 1 (i) 2 + 6(i) 2 + & 2 (x 2 (i))>0. 

Proof. As 6", u and are Lipschitz, by the Cauchy-Lipschitz theorem the system of ODEs 

(|5.1|) has a unique maximal solution. In order to prove that this solution is defined globally, it 

is enough to prove that the time derivative of this solution is uniformly bounded. This comes 

I I 

from assumption f|2.3[) giving an upper bound on b'/(l + b 2 (x 2 )) 2 and b"/(l + b 2 (x 2 )) 2 , and 
from the lower bound on to be established now. 

The crucial assumption here is the fact that b' and b do not vanish simultaneously, and more 
precisely the existence of n, f3 > such that for all x 2 € R. 

\Hx 2 )\< V ^\b'{x 2 )\>p. 

Along a trajectory of the Rossby Hamiltonian, tr is conserved. Therefore, 

• either b(x 2 ) > 77, and (£) 2 > rj 2 , 

• or 6(x 2 ) < r\. In this last case, if |£| > ??, (£) 2 > r] 2 . Else, 

77 

M ^ 77T ~ IMIoo»7 

from which we deduce that 

n 

(0b> 



T\ + \\u\\ooV 
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In any case, we have the lower bound 

(Ob > ??min ( 1, ]_ 

V l r l + IMloo??, 

which can be made uniform for initial data in C replacing \t\ by max^ \t\. □ 
5.2. Dynamics outside from the support of u. 

Using the fact that u has compact support, and simple properties of the Rossby dynamics in 
the absence of zonal flow, we can prove the following 

Proposition 5.2. Let C be some compact subset o/R 4 such that 

cn{(x 1 ,x 2 ,t 1 ,z 2 )/ei+f 2 + b 2 (x 2 ) = 0} = ®. 

Then the bicharacteristics of the Rossby Hamiltonian starting from any point ofC are bounded 
in x 2 : 

sup 1^2(01 < °°- 

Proof. • Let us start by describing the dynamics in the absence of zonal flow: £1 is then an 
invariant of the motion, so that the dynamics in (x 2 ,£ 2 ) can be decoupled. Furthermore, as 
the energy surfaces are compact 

tib'(x 2 ) 

the motion along x 2 is periodic (with infinite period for homoclinic and heteroclinic orbits). 
The motion along x\ is then determined by the equation 

dt b[X2) (01 

It is trapped if and only if the average of the right-hand side over one period is zero. Outside 
from saddle points, this quantity depends continuously on £1, so that we expect the initial 
data leading to trapped trajectories to belong to a manifold of codimension 1. This can be 
proved rigorously if b 2 has only one non degenerate critical points (see [1]). 

• Let us now turn to the influence of the zonal flow. We will first check that the only possible 
escape direction is again x±. Indeed the energy surfaces corresponding to r ^ are bounded 
in the x 2 direction : as x 2 — > ±00, 

^+s(*).{^o. 

Consider now a trajectory on the energy level tr = 0, and some point of this trajectory 
(yi, 2/2,^1,^2) such that y 2 £ Supp^u. One has 

&'(y 2 )& = 0. 

- If b'(y 2 ) = 0, then 

dx\ dx 2 d£i d£ 2 n 
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The uniqueness in Cauchy-Lipschitz theorem implies then that the trajectory is nothing else 
than a fixed point, and therefore in particular is bounded. 

- If £1 = 0, then 

dx 2 d£i d& n , dxt , ^ + b 2 (y 2 )-^ 

-dT = ir = iF = and -dT = Hm) — w b — ' 

meaning that the trajectory is a uniform straight motion along x\. In particular, it is bounded 
in the ^-direction. 

Finally, we conclude that trajectories on the energy level tr = are either trapped in the 
support Supp^u, or trivial in the ^-direction. □ 



5.3. Perspectives. As recalled in the introduction, it is generally believed that in the sit- 
uation depicted in this paper (a flow around a large macroscopic current) Rossby waves are 
trapped. However due to the 2-dimensional setting (compared to the work in [lj) the trapping 
in the x\ direction seems difficult to prove, outside some specific cases studied in the previous 
paragraph. One way to be convinced of the trapping phenomenon should be by implementing 
the dynamical system numerically. It should be pointed out however that actually in order 
to get physically relevant predictions for the oceanic eddies, one should consider 3D models, 
or at least 2D models involving the influence of stratification. The methods presented here 
seem to be robust and should be extended to such complex models, up to again the study of 
the hamiltonian system describing the Rossby dynamics. 



Appendix : a comparison result 

For the sake of completeness, we state here the result which shows the stability of the prop- 
agation under a 0{e°°) error on the propagator. This result has been implicitly used in the 
proof of the diagonalization when comparing A and T±, Tr, and in the proof of dispersion 
when comparing T± and T±. 

Proposition 5.3. Let A £ and A £ be two pseudo-differential operators such that 

• iA £ is hermitian in L 2 (R. d ), 

• A £ — A £ = 0(e°°) microlocally on C R M . 

Let tp be a solution to 

id t Cp + A £ (p = 
microlocalized in Vt, and <p be the solution to 

id t (p + A £ Lp = 

with the same initial data. Then, for all JVgN, 

sup Mt)-v(t)\\ L ^)=0{e™). 

t<e- N 
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Proof. The proof is based on a simple energy inequality and is completely straightforward. 
We have 

-J^y - ^llia(Rd) = 2(iAe<P ~ iA £ (f>\<p - <p) 
= 2{(iA £ - iA £ ){p\y - (p) 

< 2\\(A £ - i £ )(p|| L 2 (Rd) ||^ - <p\\ L 2(Bd) ■ 

This leads to 

y(t)-m\ 2 LH K* ) = o(e oo )t, 

which concludes the proof. □ 
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